Analytic model for the ballistic adsorption of polydisperse mixtures 
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We study the ballistic adsorption of a polydisperse mixture of spheres onto a line. Within a 
mean-field approximation, the problem can be analytically solved by means of a kinetic equation 
for the gap distribution. In the mean-field approach, the adsorbed substrate as approximated as 
composed by effective particles with the same size, equal to the average diameter of the spheres 
in the original mixture. The analytic solution in the case of binary mixtures agrees quantitatively 
with direct Monte Carlo simulations of the model, and qualitatively with previous simulations of a 
related model in d = 2. 

PACS number(s); 68.45.Da, 81.15.-z, 82.70.Dd, 68.10.Jy 
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I. INTRODUCTION 



The adsorption of colloidal particles onto a surface is 
a subject of considerable interest due to its many practi- 
cal applications in fields as diverse as physics, chemistry, 
biophysics, medicine, etc. |l|]. Several models have been 
proposed so far, in an attempt to understand the physical 
properties of the adsorbed phase. In the random sequen- 
tial adsorption (RSA) model [§-g|, the adsorbing parti- 
cles are located at random positions on the surface. If 
an incoming particle overlaps with a previously adsorbed 
one, it is rejected; otherwise, it becomes irreversibly ad- 
sorbed. The RSA model does not consider the transport 
of the particles, and focuses only on the excluded vol- 
ume effects. It is thus a valid approximation when the 
particles arrive at the surface purely by diffusion ||^. In 
the ballistic model (BM) R-H], when an incoming par- 
ticle fails to reach the surface directly, it is allowed to 
roll down over the previously adsorbed ones, following 
the direction of the steepest descent, until it reaches an 
equilibrium position. Particles that eventually rest on 
the surface are irreversibly adsorbed; otherwise, they are 
rejected. The BM is therefore a good approximation to 
describe adsorption in presence of strong interactions, at- 
tracting the particles towards the surface |lQ-|l^]. 

In their original formulation, the aforementioned mod- 
els, as well as their main variations, consider essentially 
the adsorption of a monodisperse suspension, in which 
the adsorbing particles have all the same size. Real-life 
suspensions, however, always possess an unavoidable de- 
gree of polydispersity. For instance, in some experimental 
situations the standard deviation of the particle size dis- 
tribution may be up to 5-10% of the mean particle size 
|^,Q. Under such conditions, the effects of polydisper- 
sity may be indeed important. 

The role of polydispersity has been studied in some de- 
tail in the RSA model. Theoretical works have dealt with 
binary mixtures of particles with greatly differing diam- 
eters ll^ , power-law size distributions ||l^,0 , or general 



continuous size distributions [|l8| . Numerical simulations, 
on the other hand, have been performed in a wider va- 
riety of conditions: binary mixtures [|l9| , p0|, u niform size 
distributions [|l9[, gaussian distributions [14 19|, power- 
law distributions ||l^, etc. 

In the context of the BM, it is noteworthy the work 
of Senger et al. iQ , where the authors describe a Monte 
Carlo model for the adsorption of a two-component mix- 
ture of hard spheres onto a plane, where the particles are 
under the simultaneous influence of diffusion and gravity. 
This is indeed a mixed model, which reproduces the stan- 
dard RSA model in the limit of small particles, and the 
BM for large particles. The results reported by Senger 
et ai, for particles large enough to be well inside the BM 
regime, are qualitatively similar to those found for the 
RSA of binary mixtures p9[ : The maximum fraction of 
surface covered by the adsorbed particles — the jamming 
limit 9oo — increases monotonically with the concentra- 
tion p of large particles, with a maximum in the limit 
p ^ 1^ {i.e., 1 — p arbitrarily small, but non-zero). For 
p = Q and p = 1, the coverage corresponding to monodis- 
perse adsorption is recovered. 

In this paper we present an analytic model for the bal- 
listic adsorption of mixtures of spherical particles with 
different diameters. The model can be solved in a mean 
field approximation, by studying the kinetics of the gap 
density function [g|. Within this approach, we are able to 
derive a generic equation for an effective gap distribution. 
To test our equation, we solve it explicitly in the simplest 
case of a binary mixture. The analytic results obtained 
for the density at jamming Ooo match the findings of di- 
rect Monte Carlo simulations of the model. Moreover, 
the qualitative behavior of 0oo predicted by our model is 
the same as the reported by Senger et al. [Q . 



II. MODEL 

Our model considers the adsorption onto a line of a 
polydisperse mixture whose degree of polydispersity is 



characterized, in general, by a continuous distribution of 
sizes p((t). The quantity p{a)da is defined as the fraction 
(bulk concentration in the infinite reservoir from which 
the particles are drawn) of spheres with diameter between 
a and a + da. We assume p to be normalized to 1. Thus, 
for a monocomponent solution of particles of size CTq, we 
have p{a) = 6(<7 — ao)- The particles arrive at the line 
at rates k{a) per unit length per unit time. Assuming 
that the adsorbed substrate interacts with the incoming 
particles only through excluded volume effects, we can 
select the appropriate units of time and set k{a) = p{a). 
Under these conditions, the problem translates into the 
sequential adsorption of particles of size a, selected with 
a probability density p{a). 
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FIG 1. Landing configurations for particles of different sizes 
(Ti and (72 > o"i. 



When an incoming particle lands over a preadsorbed 
one of exactly the same size, the adsorption rules are 
identical to the standard BM [g| . Figure [^ depicts the 
possible configurations involving particles of different di- 
ameters cTi and (72, with ui < 02- When a small particle 
rolls over a large one, the former finally falls on the sur- 
face and, after it is adsorbed, the centers of both parti- 
cles are separated a horizontal distance A = {a\ + cr2)/2; 
Fig. |l|(a) represents this case. When a large particle rolls 
over a small one, the rule adopted in our model is the 
one represented in Fig. 0(b), in which, after rolling, the 
centers of the two particles are also separated a distance 
A. 



The adoption of the rule pictured in Fig. ||(b) repre- 
sents a major simplification in the model. It could be 
possible to argue that, in a more realistic treatment, the 
final configuration involving a large particle rolling over 
a small one should be the one depicted in Fig. Kc) . The 
surfaces of the particles are tangent after adsorption in 
this case, and their centers are separated a horizontal 
distance A^ = ^a\ oi . Both rules can be easily imple- 
mented in a numerical simulation. However, the prescrip- 
tion nVc) imposes an essential asymmetry among particles 
of different size. First of all, in our model, as defined by 
rules ll|(a) and |l|(b), the final result of an adsorption event 
involving two spheres of different diameter is independent 
of the order in which the particles reach the surface. As a 
consequence, our model does not allow for "overhangs" ; 
this means that, if n((T) is the density of adsorbed par- 
ticles of size (7, then the fraction of covered surface 9 is 
simply given by = JdcrCTn(CT). This simple expression 
obviously does not hold in a model defined with the rule 

0(c). 

These "abelian" properties are eventually responsible 
for our model being analytically tractable. 



III. GENERAL MEAN-FIELD EQUATION 

The model defined in the previous section can be ana- 
lyzed by studying the density function of gaps — holes be- 
tween two consecutive adsorbed particles. Let us define 
G(a:,i)da: as the number of gaps with a length between 
X and X + dx present at time i, per unit length of sub- 
strate. The time evolution of G is obtained as a balance 
equation for the creation and destruction of gaps caused 
by a single adsorption event [gj . Given G, the fraction of 
covered surface is defined by 



d{t) = 1- xG{x,t)dx, 



(1) 



and, from here, we obtain the jamming limit as 0oo = 

liTat^ao9(t). 

In the case of the ballistic adsorption of a monodis- 
perse solution of spheres of diameter CTq, the equations 
for the density of gaps are M 



dG{x,t) 

dt 
dG{x,t) 

dt 



-{x + ao)G{x,t) + 2aoG{x + ao,t) + 2 G{y,t)dy, for x > ao] 



X'\-(To 



2(JoG{x + cTo, t) + 2 / G{y, t)dy, for x < cfo- 

J X-\-(7o 



(2) 
(3) 



The solution of (|)-(|) is 

G(x,t) = e-(^+'"°'*t2^(aoi)exp {2(1-6-'"°*)}, for x > ao; 

G(x,i) = 2 /" duu(l-hCTou)e-(^+2'"°)"F(CToM)exp{2(l-e-'"°")}, for x < ao, 

Jo 



where we have defined the auxihary function 



F{t) = cxp <^ -2 



1-e- 



(4) 



For a polydisperse mixture, the naive application of 
this approach becomes considerably more involved. Af- 
ter a moment's reflection, it is easy to realize that, in this 
case, the final configuration resulting from an adsorption 
event taking place on a given gap, depends on the sizes 
of the particles defining the boundaries of that gap. We 
should accordingly deal with a continuous set of func- 
tions Gcr',cr" {x, t), defined as the densities of gaps created 
between particles of size a' and a", for a', a" G [0, cxj] . 
An enumeration of all the possible events occurring when 
adsorbing spheres of size a at rate p{(t), would lead to 
a system of exact coupled integro-differential equations 
for the magnitudes Ga'^a", that would completely deter- 
mine the dynamics of the process. The magnitude of this 
task, especially when dealing with continuous size distri- 
butions p(cr), seems to preclude any chance for an exact 
solution. 

Fortunately, however, a great deal of insight can be 
gained by seeking for a mean-field type of solution, based 



in the following argument: When the particles are free in 
the suspension, they are distinguishable and interact dif- 
ferently with the adsorbed phase, depending on their size. 
However, once they have been adsorbed, we can assume 
that they become indistinguishable, in the sense that the 
adsorbed particles interact with the incoming particles 
as if the former were all equal, with the same average 
diameter a — J ap{a)da. In other words, we can ap- 
proximate the adsorbed phase as composed by effective 
particles of the same size a, interacting with incoming 
particles of size a. Assuming this simplification, we need 
only a single effective gap distribution G, defined by the 
gaps bounded by the adsorbed effective particles. 

We remark the important fact that the aforementioned 
mean-field approximation does not imply at all that the 
density of adsorbed particles is proportional to the bulk 
density, n{a) ex p{cr). This relation, which can be true at 
the first stages of the adsorption process, does not hold 
close to the jammed — saturated — state. This last state- 



ment is most easily seen in binary mixtures (see Sec. IV). 

The kinetic equation for the effective gap density can 
be written in the generic form 



dap{a){x + a)G{x,t)+2 da p{a) (^ + ^) G{x + cr,t) + 2 da p{a) dyG{y,t). (5) 



dG{x,t) 
di 



The origin of the different terms in (H) is the following: 
The destruction of gaps of length x is due to the land- 
ing of a particle of size a on any point of an interval of 
length X + a centered on the gap. After averaging over 
the distribution of incoming particles of size a < x, we 
obtain the first term in (|^). A gap of length x can be 
created by the impact of particles of size a on either of 
the particles of effective size a defining a gap of length 
X + a. These events, which happen at rate p{a), account 
for the second term in (|^). The last term is due to the 
averaged creation of gaps of length x by direct deposition 
of a particle of size a onto a gap of length y > x + a. We 
remark again that Eq. (g) owes its relatively simple form 
to the choice of the "abelian" rule ||(b) in the definition 
of the model. A much more complex expression would 
have been obtained with rule |l|(c) . 

Equation (^ can be expressed in a more compact way 
by integrating by parts its last term. Defining the distri- 
bution function ^(x) = L p{a)da, we obtain 



X-\-(7 



dG{x,t) 



dt 



= --^{x){x + a)G{x,t) 



da 



p{a){a + a) + 2'^{a) G{x + a,t) 



(6) 



Equation (pi) is the final expression of the mean-field 
theory for our model of polysdisperse ballistic adsorp- 



tion. As a consistency check, we consider the trivial sce- 
nario of a monodisperse suspension. In this case, by set- 
ting p{a) = 5{a — (To) and '^{x) — Q{x — CTo), where O 
is the Heaviside step function, we immediately recover 
the equations for a single-size distribution, as given by 
Eqs. d) and d). 



IV. BINARY MIXTURES 

In order to test the validity of our mean-field theory, 
we proceed now to solve explicitly Eq. (ph in the case of 
a binary mixture, composed by particles of size ai = 1, 
which adsorb onto the surface at rate 0i, and particles 
of size (72 = r > 1 , adsorbing at rate (^^ = 1 — 0i ■ As 
an aside, in this simple setting we can estimate the vari- 
ations in the jamming limit due to the adoption of rule 
|l|(b) instead of |](c). One can expect that, for small val- 
ues of r, the outcome of both models should be similar. 
Indeed, numerical simulations show that, for values of 
r < 2, the difference between prescriptions is always less 
than 1%, for all values of (t>r- 

The density function for a binary mixture has the form 
p{a) — (j)i5{a — 1) -|- (f>r5{a — r), whereas the distribution 
function is ^(x) — (j)iQ{x — l)ld{r — x) + Q{x — r), and the 
average size a = 4>i-\- r(f>r ■ By inserting these expressions 
into (H) or ffl), we obtain the following set of equations: 



dG{x,t) 
di 



-{x + a)G{x, t) + 0i(l + a)G{x + 1, t) + 0r(r + d-)G{x + r, i) 



dG{x,t) 
dt 



dG{x,t) 
dt 



+ 2(/)i/ Giy,t)dy + 2q^r G{y,t)dy, for x > r; 

= -01 (a; + a-)G'(a;, t) + 0i(l + CT)G(a; + 1, t) + 0r(?' + o)G{x + r, t) 

/•OO /'OO 

+ 2(/)i/ G{y,t)dy + 2(f>r G{y,t)dy, for 1< x < r; 

^2;4-l J x-\-r 

01 (1 + ct)G'(x + 1, t) + (/-rlr + a)G{x + r, t) 

201 / G{y,t)dy + 2cj,r G{y,t)dy, for < x < 1. 

J x-\-l J x-\-r 



(7) 



(8) 



(9) 



We observe that, for a binary mixture, one could in principle try to solve exactly the model, by determining the 
rate equations for the densities of gaps delimited by particles of size 1 and r, namely Gi^i, Gr,ri and Gi^r- However, in 
this case one would end up with a set of nine coupled equations. The simplification achieved through the mean-field 
theory is here evident. 

We consider in particular the case 1 < r < 2. To solve the kinetic equations, we seek in (0) a solution of the form 
G{x,t) — e~*^^+°"^*i?(t). With this substitution, we are led to the equation for H{t): 

dlnH 



dt 



(1 + .) + - 



(r + a) + - 



(10) 



The solution of dlQ), with the initial condition H{0) = 0, is 

H{t) = t^ exp {0i(l + a){l- e-*)} exp {0,(r + a)(l - e-''*)/r} [F{t)]^' [F(ri)]^^ (11) 

where F(t) is defined in (0). Upon substituting this result into (0), we look for a solution of this equation of the form 
G{x,t) = e~'^i(^+*)*g(a;,t). The equation determining Q is 



dQix,t) _ _^,.(.+^), dg(t) 
dt di ' 

from which Q{x,t) is obtained by direct integration, together with the initial condition Q{x,0) = 0: 

g(a;,<)=e-'^-(^+'^)*i/(t) + 0^(x + CT) / dw e-*'-(^+'^)"iJ(u). 



(12) 



(13) 
Jo 

Finally, by substituting the solutions of (0) and (||) into the appropriate range of values of x in Eq. (||) (and taking 
into account that r < 2), we can directly integrate this equation. Using Eq. (|l|), and after performing some algebraic 
manipulations, we obtain the density of adsorbed particles as a function of time: 

9{t) = / duH{u)J'i{u) + (t)r / duH{u)J'2i<Pit + <j)ru) 
Jo Jo 

+ 010r(l+O-) / duH{u) / dv J^3{(l)iV + (f)rU) + 2(f>i(f>r / du H {u) / dv J^4{(f>lV + (prU) , 











where we have introduced the auxiliary functions 



.Fi(z) 



2^ 
^ — az 



{ [20,r + (1 + 0r + 0rCT)z +((T+ 1)Z^] e '- - [2(t)r + 0r(r + &)z\ 6^"} 



^2{z) = —^ {- [2 + (a + 2)z +(a + 1)2^] e"^ + [2 + (a + 2r)z + r{a + r)z''] e""} , 



Mz) 



Z 

^ — az 



^ — az 



{- [2 + ((7 + l)z] e-' + [2 + (a + 2r - l)z + (r - l)(r + S)^^] g-''"} 



.F4(z) = -^{-[3 + (a + l)z], 



3 + (o- + 3r - 2)z + i(r - l)(3r + 2ct - 1)^^ + i(r - lf{r + d)z^ 



We can estimate the theoretical predictions of this 
mean-field solution by numerically integrating the pre- 
vious expression in the limit f — > oo. Figure shows in 
full lines the results of the integration for different values 
of r. The symbols represent data obtained from direct 
Monte Carlo simulations of the model on a line of length 
1000 with periodic boundary conditions. We observe that 
the predictions of the mean-field theory are in excellent 
agreement with the numerical simulations. 
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FIG 2. Jamming limit as a function of the concentration frac- 
tion (j>r of large particles, for different values of the diameter 
ratio r. Comparison between numerical simulations (hollow 
symbols) and the mean-field prediction (full lines). 



From Fig. y we conclude that, for <j)r < 1, the jamming 
limit is an increasing monotonic function of this variable. 
For 0J. = or 0r = 1 (only small or large particles, re- 
spectively), we recover, for any r, the prediction of the 
standard BM model, 6*^*^ ~ 0.808 ||. For small values 
of (j)r, doo grows linearly, 6oq — 0^ + a{r)(j)r, with an 
slope a{r) that increases with r. The value of the slope 
at the origin can be easily estimated by Taylor expand- 
ing the expression for 0{t). The jamming limit exhibits a 
maximum located at 0^ — *■ 1 , iu qualitative agreement 
with the findings of Senger et al. [|l^ . The actual value 
of the maximum 91^"^ (r) is an increasing function of r, 
with an apparent tendency to saturate at large r. In the 
limit 0J. -^ 1, and for r ^ 1, we can estimate the limit- 



ing value of 6'™°^(r) |19|: In this limit, the large particles 
cover first a fraction of surface 9^^^ of the line, leaving 
free a surface 1 — 9^^'^ that is afterwards covered un- 
til jamming by the small particles. The total coverage is 
therefore bounded by 9^'^'-'{r) < 9g^' -t- (1 - 9g^'')9g'^' = 
e^^{2-9g^) ~ 0.96339. Monte Carlo simulations con- 
firm this extreme, yielding the value ^JJ"^ = 0.964±0.001 
for r = 20 and Ar = 0.99. 



V. CONCLUSIONS 



To sum up, in this paper we have presented an exten- 
sion of the classical ballistic model RtM] 7 describing the 
ballistic adsorption onto a line of a polydisperse mix- 
ture of spherical particles of different sizes cr, present 
with a bulk concentration p{a). The model is solved by 
means of a mean-field-like equation for the gap density, 
which approximates the adsorbed phase as composed by 
effective particles, all with the same average diameter 
a = Jdaap{<7), interacting with incoming particles of 
variable size. To check our mean-field approximation, 
we have explicitly solved the case of a binary mixture. 
The perfect match of the theoretical solution with the 
numerical simulations confirms the validity of the mean- 
field approximation, at least for this particular case. Our 
findings agree also with numerical simulations of a related 
model in two dimension llfl]. On theoretical grounds, 
the proposed mean-field approach could be a first step 
to deal with more complex situations, as, for example, 
in higher dimensionalities, where the assumption of an 
effective layer of adsorbed particles would be more rea- 
sonable, or in the case of adsorption onto a substrate 
initially covered with impurities. 
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